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Abstract 


Reid  measurements  with  a  simple  660  nm  forward  scattering  sensor  were  analysed  by  means 
of  a  three-layered  sphere  model  of  light  scattering.  The  scattering  sensor  measures  a  value  that 
represents  the  scattering  function  multiplied  by  a  weighting  function..  The  geometry  of  the 
scattering  sensor  determines  the  weighting  function.  A  model  was  developed  that  determines  the 
weighting  function  based  on  sensor  geometry.  For  the  sensor  tested  the  weighting  function  peaks 
at  27  degrees  with  a  FWHM  of  22  degrees. 

Measurements  with  the  device  were  made  off  the  Oregon  coast  and  off  San  Diego.  The  device 
tested  displayed  segmented  linear  relationships  with  the  particulate  beam  attenuation  coefficient  at 
660  nm  and  the  particulate  total  scattering  coefficient  at  488  nm  as  a  function  of  depth.  The  slope 
of  the  correlation  with  the  particulate  beam  attenuation  coefficient  changed  by  a  factor  of  five  over  a 
depth  range  of  125  meters. 

The  volume  scattering  functions  off  Oregon  were  calculated  using  historical  data  for  particle 
size  distributions  and  particle  densities.  The  predicted  output  of  the  scattering  sensor  using  the 
calculated  scattering  functions  and  the  calculated  weighting  function  showed  that  the  segmented 
linear  relationships  can  readily  be  explained  by  means  of  a  change  in  the  shape  of  the  volume 
scattering  function  with  depth. 

.The  scattering  sensor  thus  measures  a  weighted  integral  of  the  forward  scattering  function.  Its 
output  is  dependent  on  both  the  total  particle  concentration  and  the  shape  of  the  volume  scattering 
function  in  a  predictable  way. 


Introduction 


In  order  to  rapidly  assess  particle  fields  in  the  ocean,  a  small,  inexpensive,  optical  sensor 
without  critical  alignment  problems  is  desired.  We  have  data  from  an  early  prototype  forward 
scattering  sensor  that  satisfy  these  criteria.  Briefly,  without  impinging  on.  proprietary  information, 
it  consists  of  a  light  source,  a  detector  and  a  light  block  in  between  to  keep  direct  light  from 
reaching  the  detector.  This  geometry  is  similar  to  the  old  Lamont-Daugherty  nephelometers 
(Thorndike,  1975).  This  arrangement  results  in  the  reciever  accepting  light  that  depends  on  the 
geometry  of  the  device.  In  the  case  of  the  instrument  tested  light  scattered  at  angles  between  8°  and 
60°  is  accepted  with  varying  efficiencies.  We  can  numerically  estimate  the  efficiency  as  a  function 
of  scattering  angle  by  dividing  space  up  into  volume  elements  and  summing  up  the  amount  of  light 
that  would  hit  the  detector  due  to  scattering  at  each  angle  assuming  the  volume  scattering  function 
to  be  umformly  1.0  rn'l  ster.-f  Given  a  real  volume  scattering  function  one  could  then  dot 
multiply  by  the  efficiency  vector  to  obtain  the  signal.  Effects  accounted  for  in  this  program  are  the 
relative  intensity  of  the  light  source  in  each  direction,  the  attenuation  of  light  from  the  source  to  the 
volume  element  and  on  to  the  detector,  and  the  solid  angle  subtended  by  the  detector  as  seen  from 
the  volume  element. 

Field  test  of  Newport,  Oregon 

The  first  in  situ  test  of  this  scattering  sensor  was  carried  out  on  August  14,  1990.  Casts  were 
made  at  two  stations;  one  at  5  Nautical  Miles  and  one  at  10  NM  from  Newport,  Oregon.  In 
addition  to  scattering,  we  measured  termperature  and  beam  attenuation  (660nm)  as  a  function  of 
depth.  The  first  station  (Figure  1)  at  5  NM  offshore  was  a  simple  three  layered  system  with  a 
surface  turbid  layer,  a  clearer,  mid-depth  layer  and  a  bottom  nephcloid  layer  (BNL)  with  the 
associated  gradients  between  each.  The  second  station  (Figure  2)  at  10  NM  was  more  complex 
with  an  intermediate  turbid  layer  between  the  surface  turbid  layer  and  the  clean  water. 

Relative  to  the  particulate  attenuation,  scattering  is  enhanced  at  the  bottom,  reduced  in  the 
intermediate  turbid  layer  and  further  reduced  in  the  surface  turbid  layer.  The  different  relationships 
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are  clearly  shown  in  Figure  3.  The  linear  segments  in  the  relationship  between  attenuation  and 
scattering  signal  are  suggestive  of  linear  mixing  of  water  types.  The  region  between  20  and  50 
meters  is  suggestive  of  mixing  of  three  different  water  types.  The  slopes  of  these  relationships 
vary  by  almost  a  factor  of  five.  The  same  plot  for  station  1  (Figure  4)  shows  just  two  linear 
relationships.  The  scattering  signal-attenuation  relationships  are  almost  identical  for  the  BNL's 
from  both  stations.  The  relationship  for  the  inshore  surface  waters  is  similar  to  that  for  the 
offshore  50-90  m  water.  The  magnitudes  of  these  slopes  are  a  result  of  the  mixing  process.  The 
actual  ratios  (Figure  5)  vary  only  by  a  factor  of  three. 

Anadysis 

Without  further  analysis  it  would  not  be  clear  whether  the  change  in  the  ratio  of  scattering 
signal  to  attenuation  is  due  to  a  real  change  in  b/C  or  to  a  change  in  the  shape  of  the  volume 
scattering  function  resulting  in  the  sensor  measuring  a  different  portion  of  total  scattering.  Also,  a 
real  increase  in  b/C  can  be  due  to  either  an  increase  in  pigment  absorption  per  cell  or  a  change  in 
size  distribution  or  tnore  likely  both. 

The  easiest  question  to  answer  is  whether  a  change  in  the  volume  scattering  function  could 
acount  for  the  observed  behavior.  We  are  very  familiar  with  the  partcle  size  distributions  (PSD)  in 
the  coastal  upwelling  region  (Kitchen  et  al.,  1975,  1978;  Pak  ct  al.,  1980).  In  surface  waters  PSD 
slopes  are  low  in  nutrient  rich  waters  and  high  in  relatively  nutrient  poor  water.  The  high  values  of 
C  and  the  fact  that  the  maximum  is  right  at  the  surface  indicates  that  we  are  dealing  with  nutrient 
rich  waters  and  thus  low  PSD  slopes.  The  intermediate  waters  just  below  the  surface  layer  are 
known  to  have  average  slopes  and  the  BNL  is  known  to  have  high  slopes.  We  also  know  that  the 
apparent  density  of  the  particles  (dry  weight/wet  volume)  increases  with  depth  from  1.0  g/cc  or 
less  at  the  surface  to  2.0  or  more  in  the  BNL  (Peterson,  1977).  From  this  we  may  infer  an 
increase  in  the  index  of  refraction  of  the  particles.  The  fact  that  the  near-surface  light  attenuation 
length  is  less  than  half  a  meter  also  suggests  that  the  chlorophyll  content  of  the  cells  should  be 
nearly  maximum. 
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Thus  we  will  model  the  surface  layer  with  a  Junge  type  size  distribution  dN  =  NoD'^dD  with 
the  "slope"  (s)  set  to  3.0.  To  model  the  diatoms  which  are  dominant  in  such  a  system,  we  use  a 
three-layered  sphere  (Mueller,  1974;  Kitchen  and  Zianeveld,  1992)  with  a  high-index  of  refraction 
outer  shell,  a  chloroplast  layer  also  with  a  relatively  high  index  which-varies  according  to  the 
anomalous  dispersion  relations  and  a  low-index  inner  core.  For  the  intermediate  waters  we  will 
use  the  same  structure  but  with  a  slope  S  =  4.0.  For  the  BNL  we  will  use  a  homogeneous  sphere 
with  a  uniform  index  of  1.09  -  O.OOSi  and  a  slope  S  =  5.0.  The  resulting  volume  scattering 
functions  normalized  to  particulate  attenuation  are  shown  in  Figure  6. 

The  BNL  has  a  much  larger  portion  of  its  scattering  in  the  8-60°  region,  so  that  the  scattering 
sensor  will  receive  more  light  for  the  same  total  scattering  coefficient.  The  mid-depth  model  also 
has  significantly  more  scattering  in  this  region  than  the  near-surface  model.  Halving  the  pigment 
content  of  the  near-surface  model  increased  the  portion  of  scattering  in  this  range  only  slightly,  A 
sensitivity  function  for  this  instrument  or  a  very  similar  one  (we  no  longer  have  access  to  design 
data)  is  shown  along  with  a  linear  blowup  of  the  volume  scattering  functions  in  Figure  7.  The 
numerical  results  of  meshing  the  sensitivity  function  with  the  volume  scattering  functions  are  given 
in  Table  I. 


Table  I.  Predicted  ratios  for  the  various  models  of  suspended  particulates  in  the  coastal 
upwelling  region.  BNL  refers  to  the  model  for  the  bottom  nepheloid  layer;  INT  to  the  model  for 
the  intermediate  layer  and  SUR  to  the  surface  layer.  HLF  is  the  same  model  as  the  surface  layer 
with  half  the  pigment  concentration.  SL#  indicates  a  slope  of  -#  for  the  differential  particle  size 

distribution. 


BNL  SL5 

INT  SL4 

HLF  SL3 

SUR  SL3 

0.362 

0.181 

0.135 

0.127 

0.313 

0.163 

0.124 

0.1 10 

b/c 

0.866 

0.903 

0.920 

0.863 

The  model  gives  almost  identical  results  to  the  observations.  That  is,  the  ratio  of  scattering 
signal  to  total  attenuation  is  almost  three  times  as  high  in  the  BNL  as  at  the  surface  and  twice  as 
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high  at  mid-depth  as  in  the  surface  layer.  It  should  be  noted  that  the  ratio  of  b/c  varies  only  by  a 
few  percent  even  though  the  size  distribution  and  pigment  content  were  varied  a  great  deal.  We 
can  say  that  changes  in  the  particle  size  distribution  and  index  of  refraction  are  quite  likely  to  give 
the  results  shown  for  the  field  data.  That  is  not  to  say  we  can  dismiss  the  other  factors. 
Knowledge  of  the  size  distnbution  below  2  diameter  is  poorly  known.  At  submicron  sizes, 
changes  in  the  pigment  content  can  increase  the  beam  attenuation  with  little  influence  on  the 
forward  scattering  and  thus  change  b/C  without  changing  the  forward  portion  of  the  volume 
scattering  function.  Also  our  models  may  not  realistically  proportion  the  chlorophyll  amongst  the 
various  size  classes. 

Field  data  from  off  San  Diego. 

This  scattering  sensor,  or  a  similar  one,  was  deployed  off  San  Diego  in  August  1991.  Also 
deployed  were  a  fluorometer,  standard  transmissometer  (660  nm),  and  prototype  absorption  and 
attenuation  meters  (488nm)  (Zaneveld  et  al.,  1992).  Station  7  was  the  best  example  of  the  complex 
structure  (Figure  8)  we  found.  There  are  three  distinct  particle  maxima  in  the  top  100  meters. 
Only  the  middle  one  is  associated  with  a  significant  fluorescence  peak.  However,  the  scattering 
sensor  follows  the  beam  attenuation  for  both  the  middle  and  deep  maxima,  but  is  greatly  reduced  in 
the  shallowest  maximum.  The  absorption  (488nm)  follows  the  fluorescence  very  well  (Figure  9) 
except  thar  it  gives  a  slight  indication  of  the  third  maximum  while  the  fluorescence  does  not.  The 
absorption  is  also  slightly  reduced  at  the  surface  relative  to  the  deepest  waters,  while  the 
fluorescence  is  slightly  elevated  relative  to  the  deepest  waters. 

Since  we  have  both  absorption  and  attenuation  coefficients  at  488  nm  we  can  derive  the  total 
scattering  coefficient  by  subtraction  and  compare  that  to  the  measurement  from  the  660  nm 
scattering  sensor  (Figure  10).  b(488)  shows  all  three  maxima  strongly  similar  to  c(660)  while  the 
scattering  sensor  measurement  at  660  nm  is  only  half  as  strong  in  the  shallow  maximum.  Scatter 
plots  of  Cp(488)  and  bp(488)  versus  Cp(660)  (Figure  11)  show  one  linear  relationship  for  the 
entire  cast  for  bp(488)  vs  Cp(660)  and  just  a  slight  reduction  of  Cp(488)  vs  Cp(660)  in  the  shallow 
maximum.  The  approximately  10%  reduction  in  bp(488)  can  easily  be  attributed  to  changes  in  tlic 
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absorption  spectrum.  However,  a  plot  of  the  scattering  sensor  signal  (660  nm)  versus  the  derived 
total  scattering  coefficient  at  488nm  (Figure  12)  shows  a  separate  linear  relationship  for  the  shallow 
maximum  compared  to  the  rest  of  the  cast.  The  scattering  sensor  again  showing  half  the  signal  in 
the  shallow  maximum  as  in  the  rest  of  the  cast.  Since  the  absorption  and  fluorescence  is  low  here 
it  would  be  hard  to  attribute  this  to  anomalous  dispersion  or  refraction  effects.  Thus,  the  data  from 
this  cruise  gives  even  stronger  evidence  that  the  anomalies  in  the  scattering  sensor  are  due  to 
changes  in  the  volume  scattering  function  and  are  not  reflections  of  the  behavior  of  particle 
concentration  or  total  scattering  coefficients. 

It  would  have  been  interesting  to  simultaneously  deploy  scattering  sensors  weighted  more 
towards  smaller  angles,  which  was  in  fact  attempted.  However,  the  npwer  sensors  failed  each  time 
and  thus  we  are  left  with  only  data  from  the  older  design. 

Conclusions 

1)  The  relationship  of  the  output  of  the  scattering  sensor  can  be  modeled  by  means  of  an 
analysis  of  its  geometry  and  a  knowledge  of  the  nature  of  the  particles  being  measured. 

2)  A  scattering  sensor  measuring  at  the  middle  forward  angles  cannot  be  expected  to  yield 
more  than  a  qualitative  picture  of  the  particle  concentration  or  total  scattering  coefficient  profiles 
due-to  variations  in  the  shape  of  the  volume  scattering  function. 

3)  Scattering  sensors  may  provide  extremely  interesting  information  if  the  appropriate 
ancillary  measurements  are  taken  to  make  an  interpretation  of  the  signal  possible. 
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re  1.  Profiles  of  the  beam  attenuation  coefficient  at  660  nm  and  the  output  of  the 

scattering  sensor  (V)  as  a  function  of  depth  taken  on  August  14,  1990  5  NM  off  Newport. 


Figure  2.  Profiles  of  the  beam  attenuation  coefficient  at  660  nm  (m'^)  and  the  output  of  the 
scattering  sensor  (V)  as  a  function  of  depth  taken  on  August  14,  1990  10  NM  off  Newport, 


Scattering  signal  (volts) 


Particulate  attenuation  (1/m) 

Figure  3.  Scatter  plots  and  regressions  for  various  depth  intervals  of  the  particulate  beam 
attenuation  coefficient  at  660  nm  (m"^)  and  the  output  of  the  scattering  sensor  (V)  taken  on 
August  14,  1990  10  NM  off  Newport,  Oregon. 
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Particulate  attenuation  (1/m) 

Figure  4.  Scatter  plots  and  regressions  for  various  depth  intervals  of  the  particulate  beam 
attenuation  coefficient  at  660  nm  (m“^)  and  the  output  of  the  scattering  sensor  (V)  taken  on 
August  14,  1990  5  NM  off  Newport,  Oregon. 
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Scattering  signal/particulate  attenuation  (volts  m) 

Figure  5.  Ratios  of  the  output  of  the  scattering  sensor  (V)  and  the  particulate  beam  attenuation 
coefficient  at  660  nm  (m'^)  for  tJie  two  stations  shown  in  fiourcs  1  and  2. 


Scattering  angle  (degrees) 


Figure  6.  Volume  scattering  functions  calculated  with  a  three  layered  sphere  model  for  the  Bottom 
Nepheloid  Layer  (BNL),  the  intermediate  water  (Full  pigments  s=3),  the  surface  water  (Full 
pigments  s=4)  and  the  surface  water  with  half  of  the  pigments  (  Half  pigments,  s=:3).  For 


details  see  text. 


Volume  scattering  funct./cp  ( 1  /ster.) 


Scattering  angle  (degrees) 

Figure  7.  The  calculated  relative  sensitivity  of  the  scattering  sensor  and  the  volume  scattering 
functions  shown  on  figure  6. 
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Depth  (m) 


Figure  8.  Profiles  of  stimulated  fluorescence  (relative  units),  the  particulate  beam  attenuation 
coefficient  at  660  nm  (m'l)  and  the  output  of  the  scattering  sensor  (V)  as  a  function  of  depth 
taken  during  August ,  1990  20  NM  Southwest  of  San  Diego,  Cal. 
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Figure  9.  Profiles  of  stimulated  fluorescence  (relative  units),  and  the  particulate  absorption 
coefficient  at  488  nm  (m‘^)  as  a  function  of  depth  taken  during  August  ,  1990  20  NM 
Southwest  of  San  Diego,  Cal. 


5 


Figure  10.  Profiles  of  the  particulate  total  scattering  coefficient  at  488  nm  (m'^)  and  the  output  of 
the  scattering  sensor  (V)  as  a  function  of  depth  taken  during  August ,  1990  20  NM  Southwest 


of  San  Diego,  Cal. 
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Figure  1 1.  Upper  panel.  Scatter  plot  of  the  particulate  beam  attenuation  coefficient  at  488  nm  as  a 
function  of  the  particulate  beam  attenuation  coefficient  at  660  nm  for  data  taken  during  August, 
1990  20  NM  Southwest  of  San  Diego,  Cal.  Lower  panel.  Scatter  plot  of  the  particulate  total 
scattering  coefficient  at  488  nm  as  a  function  of  the  particulate  beam  attenuation  coefficient  at 
660  nm  for  data  taken  during  August ,  1990  20  NM  Southwest  of  San  Diego,  Cal. 
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Figure  12.  Scatter  plot  of  the  output  of  the  scattering  sensor  at  660  nm  as  a  function  of  the 
particulate  total  scattering  coefficient  at  488  nm  for  data  taken  during  August,  1990  20  NM 
Southwest  of  San  Diego,  Cal. 
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of  light  scattering.  The  scattering  sensor  measures  a  value  that  represents  the  scattering  function  multiplied  by  a  weighting 
function.  The  geometry  of  the  scattering  sensor  determines  the  weighting  function.  A  model  was  developed  that  determines  the 
weighting  function  based  on  sensor  geometry.  For  the  sensor  tested  the  weighting  function  peaks  at  27  degrees  with  a  FWHM 
of  22  degrees. 

Measurements  with  the  device  were  made  off  the  Oregon  coast  and  off  San  Diego.  The  device  tested  displayed  segmented 
linear  relationships  with  the  particulate  beam  attenuation  coefficient  at  660  nm  and  the  particulate  total  scattering  coefficient  at 
488  nm  as  a  function  of  depth.  The  slope  of  the  correlation  with  the  particulate  beam  attenuation  coefficient  changed  by  a  factor 
of  five  over  a  depth  range  of  125  meters. 

The  volume  scattering  functions  off  Oregon  were  calculated  using  historical  data  for  particle  size  distributions  and  particle 
densities.  The  predicted  output  of  the  scattering  sensor  using  the  calculated  scattering  functions  and  the  calculated  weighting 
function  showed  that  the  segmented  linear  relationships  can  readily  be  explained  by  means  of  a  change  in  the  shape  of  the 
volume  scattering  function  with  depth. 

The  scattering  sensor  thus  measures  a  weighted  integral  of  the  fonward  scattering  function.  Its  output  is  dependent  on  both 
the  total  particle  concentration  and  the  shape  of  the  volume  scattering  function  in  a  predictable  way. 
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